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1. INTRODUCTION 



A new aspect of the problem of hydrogen atom in constant electric {£) and magnetic {7i) external 
fields was observed recently. It was found [1], that the asymptotic of the perturbation series in 
powers of £ considerably changes at some Ti. values. This change is related with complex solutions 
of classical equations of motion, previously not taken into account. A look at this phenomenon 
from the expansion in powers of Ti point of view could be useful. Our work is aimed just in this 
direction. An effective method for building the perturbation series is offered and the asymptotic of 
this series for the Zeeman effect is discussed here. 

The moment method for high orders of the perturbation theory evaluation, possessing additional 
possibilities as compared with the other known recurrent methods, was introduced in Ref. [2]. It 
was described as an expedient for dimensional expansion investigation — in a problem leading to 
effective isotropic anharmonic oscillator. Then it was applied to dimensional expansion for three 
body problem [3], where the effective anharmonic oscillator is anisotropic. 

We turn our attention to the fact, that the perturbation of the hydrogen atom by a potential of 
polynomial form is also convenient to investigate by the moment method. The constant homoge- 
neous electric and magnetic fields are referred just to this kind of perturbations. The advantages of 
the moment method are clearly seen in the Zeeman effect problem. High orders of the perturbation 
theory (PT) using different approaches were studied in Refs. [4-6]. The method, based on the group 
theory was applied. 36 coefficients of the hydrogen's ground state energy expansion in powers of 
Ti^ were published in Ref. [6]. For excited states Zeeman's PT coefficients not higher than of 
the third order in Ti^ arc given in Refs. [7-10]. The variables in the Schrodingcr equation for the 
Zeeman effect can not be separated and it makes computation of higher PT orders more difficult. 
The moment method does not require variables separation. Besides, as it will be shown here, this 
method can be applied to degenerate states. Obviously this possibility is essential for most of the 
hydrogen atom states. 

Apparently the logarithmic perturbation theory (LPT) was most frequently used for recurrent 
evaluation of PT coefficients up to now [11-15]. Owing to its simple algebraic structure, LPT 
allowed to compute the highest orders corrections for some cases. Thus, for the Stark shift of 
the hydrogen ground state 160 orders of PT were obtained this way [16]. But LPT is not free 
from restrictions. Even one node of the wave function leads to considerable more complicated 
computation procedure [17]. Much more complicated LPT looks in problems, where variables can 
not be separated. Without variables separation only a few initial PT orders were computed with 
the help of LPT for the hydrogen atom in electric and magnetic fields [9,10,14]. 

It is also worth mentioning here the old PT version, based on the generalized virial and 
Hellmann-Feynman theorems, which anticipated the modern moment method. In Refs. [18,19] 
recurrence relations were written for problems with spherical symmetry, allowing to compute en- 
ergy eigenvalues corrections for states with nodes in as simple way as for nodeless states. But the 
moments introduced in these papers were diagonal ones and that is why the region of applicability of 
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the method is restricted. For instance, in this formahsm it is impossible to obtain the corresponding 
wave function. 

The recurrence relations of the Ader moment method will be obtained below. Two examples 
illustrate their applications: one is referred to non-degenerate state and the other to the case 
of degeneracy. Then the asymptotic of the numerically obtained energy expansion coefficients is 
considered and the results of summation of perturbation series are given. It will be shown on example 
of the ground state, how the wave function of the perturbed hydrogen atom can be obtained within 
the moment method. 



2. MAIN RECURRENCE RELATION 



Consider the state of the hydrogen atom I'^o) = with principial quantum number n, angular 
momentum I and zero projection of the angular momentum m — 0, perturbed by strong magnetic 
field H — Hz- Let us write expansions of the energy eigenvalue and the wave function of this state 
in the form 

oo oo 

E^J2Ek7'\ ^(r) = E^^(^V\ (1) 

fe=0 k=0 

7 = n^n/Ho , Ho = e^m^c/h^ = 2.35 x lO'^G . 

Here V'fc is the correction of the order k to the Coulomb wave function, which satisfies the inhomo- 
geneous equation 

k 

{Ho - Eo)iPk = -Hi^k-i + E ^j'l'k-j , (2) 



where 



(We use atomic units.) To change the differential equation by an algebraic one we introduce the 
moments of the order k 

P^, = {iPoK-'-z'-m , (3) 

where 

I -00) — Ce"'"/", a and u are integer . 

In this definition [■^o) contains only exponential factor of the unperturbed wave function, bearing 
its scale and having no nodes. (Note, that the common normalization factor of all moments can be 
chosen arbitrary.) Just as it was done in Ref. [2], multiply Eq. (2) from the left by {t/jolr'^"'^ and 
use the possibility for the hamiltonian to act to the left, on explicitly known functions. This way 
the recurrence relation for moments of the order k results: 



2 ^a-2,u H 2 ^<T-2,u-2 ^a-l,u + ^<ju^k - ricju , [V 
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where 
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>fe-l 
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The right-hand side of Eq. (4) contains moments only of preceding orders. The coefficient, 
which we call the hyper-susceptibility of the order /c, can also be expressed through the preceding 
orders moments. This expression follows from (4) and will be written down below, separately for 
each of the cases under investigation. 

Quite similar it is possible to consider another perturbation of the same state, if this perturbation 
has the form of a polynomial in r and z. It is enough to change only the right-hand side of Eq. 
(4) for this aim. The expression in brackets in the right-hand side, representing the magnetic field 
contribution, should be replaced by another function of the preceding orders moments, created by 
the new perturbation. 

The succession of the based on relation (4) computations becomes more lucid if one represents 
on a plot the lattice of points with integer coordinates of columns a and rows v. The indices of 
moments, necessary to compute energy and wave function corrections, are placed on this lattice in 
the sector a > f — 1, i/>0. In the general case Eq. (4) relates moments of the order k, indices 
of which are located in the vertices of a rectangular triangle (see examples in Fig. 1). When one 
vertex appears to be on the line a — v — 1 and another one is outside the above indicated sector, 
the triangle (example A) transforms into segment (example B), and Eq. (4) turns into relation 
between two moments from different rows. Equation (4) relates moments in pairs also along each of 
the lines with v — Q and v — 1 (examples C and D). At A; = 0, as the direct integration indicates, 
among the moments {'^Q\r'^~^ z^\n, I) are equal to zero those, for which I — l<a<n — 2. If the 
perturbation is even, as in the case of the Zeeman effect, all corrections to the wave function have 
the same parity. Therefore, in all PT orders the moments with odd sum u -\- 1 vanish. 



Magnetic field docs not mix states with different parities, therefore besides the ground state, 
2s-, 2p- and 3p-states should also be considered as non-degenerate. Let us show, how the moment 
method works in the last case. Only the moments with odd u values can be different from zero. In 
initial PT order 



For all following orders an additional requirement is introduced: corrections to the function j'^o) 
should be orthogonal to the function j'^o) itself, i.e. (V'olV'it) — ^o,k- This condition is the routine 
element of the Brillouin-Wigner perturbation theory (see, e.g. [20]). It was used in Ref. [2]. In the 
case of 3p-state it takes the form of the following additional relation between the moments 



3. ISOLATED STATES 




ae is integer . 



(5) 



pfe 1 
6 




(6) 
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To obtain an expression for E'^, we substitute into recurrence relation (4) first u — 1, a — 1, then 
1/ — 1, a — 2. The solution of the obtained system of two linear equations is 

E, = RiT' - ^-R'^^' = liPt' - - ^{Pt' - Pt') , (7) 

Po\ = 5<-' - ■ (8) 

The sum, containing hyper-susceptibilities of preceding orders, dropped out of the final expression 
for Ek owing to orthogonality condition (6). 

Equations (4), (6) and (7) form the closed system of recurrence relations. In each order k > 1 
the sequence of computations is arranged as follows. First, the coefficient E^ is evaluated with the 
help of (7). On the next step u = 1 and a = 3 are substituted into recurrence relation (4). In 
this case together with orthogonality condition (6) it forms a system of equations from which initial 
elements of the row of moments with u — 1 are obtained: 

Pii = liRu' - PsiEk) = Ip,\ . 

Successively increasing a by one, it is not difficult to come to the necessary boundary moment of 
this row. Substituting then in Eq. {4) u — 3 and cr = 3, we get the initial moment of the next row 

P33 = 3(3Pi\ - i?33"^ + P^^Ek) 

and so on. 

The boundary moments, i.e. the moments with the highest for the given order k indices a and 
ly values are determined by the following conditions. To compute the hyper-susceptibility of high 
order K the following moments are required: Pj^ from the region u — 1 < a < 3K, < u < 2K, 
then P^^ from the region u — 1 < a < 3K — 3, < v < 2K — 2 and so on. The computation of Ek 
coefficients for the other isolated states goes a bit more simple. 



4. DEGENERATE STATES IN THE MOMENT METHOD 



We consider as an example a pair of splitted by the magnetic field states \3s) and \3d). Taking 
into account degeneracy and using the functions 

|'0o) — cosq;|3s) -|-sinQ;|3d) , |'0o) ~ e~'^^^ , 

it is not difficult to obtain the zero order moments 

{a + 2)1 fSy ( ae(a + 3)(a + 4) \ 
^ 54(2a. + l) [2) i^^^ + 2^3 ^ J ' 

where 

C^V2tga . 
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The moments in odd rows are equal to zero in all PT orders, -P^^2ae+i — 0- orthogonality 
condition (V'olV'fe) — ^o,k is equivalent to the following relation between the moments: 

12P^o - ^(4 + OP20 + iP22 = I8P00 , k>l. (10) 

There are two independent ways for expressing hyper-susceptibility through preceding orders 
moments. 

a) Equation (4) at 1/ = and (7 = determines the moment P^i qi which is used on the next step 
in the system of linear equations. Note that Ek coefficient drops out of Eq. (4) as a consequence 
of Pqo = equahty. The system of equations, containing E^ arises if one put in Eq. (4) first v — 

cr — 1, and then 1/ — 0, a — 2. Its solution is 

Pqo = -9-^00"^ + 6-^10"^ ~ S'^^o'^ ■ (12) 

The obtained as a by-product moment Pqq is substituted into the right-hand side of orthogonahty 
condition (10). 

b) Substituting into Eq. [A] v — 2 and a — 2 together with the obtained Pqo value, we get the 
second independent expression for E}^ 

= ^(9i?oV' - ^Rio' + + R22') ■ m 

Equations (11a) and (lib) result in single magnetic susceptibility Ei value at two ^ values: 

It is a natural result, which in Rayleigh-Schrodinger PT follows from the secular equation. In 
what follows the notation 3s is kept for the state with small admixture of d-wave and ^ — ^1 — 
(3-\/4T— 13)/10, and 3d denotes the orthogonal to this state combination of s- and d- waves. In the 
following approximations the unambiguity condition, applied to hyper-susceptibility of {k -\- l)th 
order, fi'^+i ~ -^i+i equivalent to relation between seven unknown moments of the order k: 

mP2, - p^2) - 6(^3^0 - ^32)] - ^(1 - '^i)P!o + ^(1 + i)p',2 - PL = 

8[e(95o'o - e^i'o) - ^(1 - ^OS'20 + S'22] ■ (14) 

k 

Here S^^ = J2 EjP^^^^ and ^ = ^12- One more constraint on the moments is orthogonal- 

j=i 

ity condition (10). To obtain a closed system, Eqs. (10) and (14) should be supplemented by 
seven equations, following from recurrence relation (4). The set of unknown variables includes 
P^Q, P2Q, Pio, Pio, P22: P32: P42: Pii: Pu- 1^ is cuough to determine only two moments Pi^ and P2^2 
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from the system of nine equations. Then, with the help of already known moments and relation (4) 
it is not difficult to compute successively all necessary moments of the given order, passing line by 
line the lattice of indices, like in the case of 3p-state. 

Quite similar, it is possible to accomplish the computation for a state with arbitrary n value 
and zero projection of angular momentum. The unperturbed wave function has definite parity and 
contains g degenerate in energy terms. Therefore g — 1 independent mixing parameters explicitly 
enter the zero order moments. There are two groups of moments in every order k > 1. Recurrence 
relation (4) connects the moments P^^ with c < n — 2 and separately the moments of the same 
order but with a > n — 2. Moments from different groups are connected by recurrence relation only 
through the moments of preceding orders. 

To obtain all independent expressions for E^, one should substitute into recurrence relation (4) 
the successively increasing u values of the given parity and corresponding set of a values: 

0<u <a <n-l . 

At every u value the unambiguously solvable system of equations is obtained. Its solution contains 
independent expression for Ek and a set of moments of the order k, to be substituted in analogous 
system at the next u value. Thus, g independent expressions for Ek in terms of the preceding 
orders moments result. The unambiguity condition of Ei determines g — 1 mixing coefficients. 
The unambiguity condition of E^+i at k > 1 allows to express the moments P^^, from the domain 
a > n — 2 through preceding orders moments. The equations expressing the unambiguity of E^+i 
are supplemented by orthogonality condition {ipoli^k) — and by necessary number of equations, 
obtained from recurrence relation (4), to close the system. 

5. RESULTS 

Energy eigenvalues 

For all levels with n < 3 and m — with the help of the moment method we have obtained 
Zeeman's hyper-susceptibilities Ek up to 75th order, see Table 1. All computations were carried 
out with 32 decimal digits. Complete agreement is observed with the results of Ref. [8] and Ref. 
[7], containing first five coefficients E^ for the ground state and three initial coefficients for both 
2s- and 2p-states in the form of rational fractions. In Ref. [9] a difference was detected between 
the obtained in this work expression for E2 coefficient and its value at Z = 1 in Ref. [7]. This 
deviation is confirmed. As it follows from expression of Ref. [9], = -45.556, but our result is 
E2^^^ — —42 in agreement with [7]. 

Energy eigenvalues £'(7) of six states, obtained by corresponding power series (PS) summation 
with the help of Fade approximants [L/L](7^) and [L/L — 1](7^) are shown in Figs. 2a-2c. These 
figures represent also the region of convergence of Fade approximants. Without expansion in Ti,^, 
by means of the splines method, which is one of modifications of the variational method, energy 
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eigenvalues of states under consideration were computed for some H values in Ref. [21]. They are 
also indicated in Figs. 2a-2c. Reference [21] has the best precision among all non-perturbative 
calculations and contains a comparison of a large number of previous computations. The precision 
of PS summation with the help of Pade approximants is high enough. At 7 ~ 1 for 2s- and 3p-states 
three or four stable digits of energy eigenvalue are established, and at least two decimal digits for 
energies of the other states arc obtained. Results of PT series summation together with some results 
of Ref. [21] are represented in Table 2. This table also illustrates convergence of Pade approximants 
we used. The convergence is sharply increased with 7 decreasing, and at 7 < 0.3 the precision of 
PS sum exceeds the precision of variational calculations [21]. 

The PT coefficients approach, as the order increases, to the asymptotic, the leading term of 
which for the Zeeman effect is [4,5]: 



= i-ir'-^, i^T + 2n - 1 + ^) . (15) 



For the levels under consideration 

Du = 32 , = 128 , D2p = 64 , 





13 \^/^ 



where 

+ _ ^^^^^ 

This result was obtained by the method introduced by Bender and Wu [22]. First, the penetreabihty 
of the barrier at imaginary magnetic field value was computed in quasiclassical approximation, then 
the dispersion relation in was applied. Stricktly speaking, according to the conditions of Ref. 

[5], expression (15) is not referred to the case of 3p-state. But it is reasonable to suppose that 
formula (15) describes all six discussed here states with coefficient D^p fitted by comparing and 
Ek of sufficiently high orders. The result is D^p — 2^^/3^. The approach of exact coefficients E^ to 
asymptotic (15) is illustrated by Fig. 3. 

In Refs. [4,5] the corrections to the asymptotic E^ were obtained for some states, among which 
2s-, 2p- and 3p-states were absent. Writing the corrections as 

S^ = '=» + S+(2SF"'-' 

it is not difficult, following the method of Ref. [22], to obtain q coefficients for all missing in [4,5] 
states, see Table 3. As the number of q coefficients included in Eq. (16) is increased, the precision 
of their determination increases as well. 

For 3p-state this stability of the power correction coefficients confirms that the leading term Ek 
of the asymptotic is determined correctly. Note, that due to dispersion relation in H'^ [5] coefficients 
Ci are related with the corrections to quasiclassical approximation for the barrier penetreabihty at 
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T-C^ < 0. The straightforward computation of quasiclassical corrections is a comphcated enough 
problem. 

Wave functions 

By analogy with the anharmonic oscillator [2] case, the correction to the Coulomb wave 
function has the form of a polynomial in r and cos^, multiplied by |'0o)- The perturbation is 
a polynomial and the operator in the left-hand side of Eq. (2) does not change the suggested 
structure of \'^k)- For the ground state 

m = (E E o!ffr^cos'= em ^ (E E ^fr^-'^^'^m ■ (17) 

i=0 i=2j j=0 i=2j 

A remark about the origin of the summation boundaries should be done. The 6 dependence is 
introduced only by the expression r^cos^O, contained in Hi, therefore minimal power of radius in 
the internal sum of expression (17) coincides with the power of cos 6, and the highest power of cos^ 6 
coincides with PT order. It is possible to check with the help of Eq. (2) that the highest power 
of radius in {ipk) is bigger by three than that in A system of linear equations determining 

coefficients a\j^ follows from (17) and looks like 

k 3fc 

E E ^i+aJ+P^^ij = > (18) 

j=Q i=2j 

2p<a<3k , Q<f3<k. 

In place of the indicated a and f3 values one can choose another their set, resulting in (2k + l){k + l) 
independent equations. This possibility is useful to check the precision of computations. Wc checked 
also the orthogonality of the obtained corrections (17) to the function In our computation the 
orthogonality condition was preserved with a reasonable precision up to 18th PT order. Corrections 
to the wave functions up to second order are given in Rcf. [7]. By comparing this work and our one 
the single point of deviation was observed. The sign at the r^7^ term in Ref. [7] is erroneous and 
as a result the first order correction is not orthogonal to unperturbed wave function in work [7]. 

Within the moment method we have computed |-(/'(0)p values by summing with the help of Pade 
approximants the PT series for the normalization factor and for the wave function itself. Results 
are represented in Fig. 4. One can see that 18 PT orders allow to advance up to 7 ~ 0.4. Besides, 
there is an agreement with Ref. [23] , where the wave function of the ground state was computed 
by a different method. 

6. CONCLUSIONS 

A good agreement between energy eigenvalues obtained by PT series summation and the cor- 
responding results of independent variational calculations provides a twofold verification. First, it 
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is confirmed that the computed PT coefficients are correct. Then, the absence of non-perturbative 
contributions to energy eigenvalues is indicated (the summation method is chosen correctly). 

Our investigation demonstrates the high efficiency of the Ader moment method. Let us stress 
once more the properties of the method which allowed to obtain the high PT orders for the Zeeman 
effect. These properties could be useful for solving some other problems of the perturbation theory. 

- The simple enough recurrence relations are obtained without variables separation in the initial 
equation. With their aid the corrections computation is reduced to purely algebraic procedure. 

- The nodes of the wave function are not shown explicitly and do not complicate the computa- 
tions. As a consequence the excited states can be considered. 

- The level degeneracy, as was demonstrated above, does not restrict application of the mo- 
ment method, but it leads to more complicated sequence of operations when solving the recurrence 
relations. 

- Besides, the set of moments used for hyper-susceptibilities determination contains information 
about the wave function. Corrections to Coulomb wave function, just as corrections to energy 
eigenvalues, follow from purely algebraic procedure. 

As it was already mentioned, application of the moment method to the hydrogen atom is not 
restricted by the Zeeman effect case. Any perturbation of the kind V{r) = J2 bijr'-'^z^ leads to 
recurrence relation (4) with an obvious simple modification of the right-hand side. For example, 
this way it is possible to take into account the joint influence of external electric and magnetic flelds, 
homogeneous or consisting of finite number of multipoles. 
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Table 1. Hyper-susceptibilities of hydrogen atom in magnetic field. 



Is — state 


2s — state 


u 


Ek 




1 


+2.50000000000000000000000000000 X 


1 n— 1 
10 ^ 


+O.0UUUUUUUUUUUUUUUUUUUUUUUUUUUU X iU 


2 


—2.76041666666666666666666666667 x 


1 A — 1 
10 


— 1.09000000000000000000000000000 X 10 


o 


+1.21115451388888888888888888889 x 


10 


1 O OJ::.nO A A A A A A A A A A A A A A A A A A A A A A A A A 1 a4 

+z.zo0o4444444444444444444444444 X 10 


A 

4 


—9.75540590639467592592592592593 x 


10 


— o.oiDbo44zyDzyDzyDzyDzyDzyDzyDo x lo 


5 


+1.17863024612238377700617283951 x 


1 a2 
10 


+l.oolDoOyzSlOz / lD04yooz/ 1004900 X 10 


D 


—1.95927276058352435076678240741 x 


1 a3 
10 


O O A/I O A1 TOi^ 1 OA O AOi^ yl 1 nTC O AOi^ /I 1 All 

— o.z04oyi /oDolz0ooo0oD4iy /oo0oD4 X 10 


7 


+4.27486169952196866486628589997 x 


10 


1 A QQnc;ni^'70 C^TOQi^AAA/l AOC;OOQQ C^i^nTC^ 1 a14 

+4.009090 /00/Z5DUUy4Uzozz550Dy/0 X iU 


o 

o 


—1.18693528256085740621952195187 x 


10 


O O 1 O AT/^TAOAAC:/^ Ai^OOAAyl yl yl O 1 OO AA x^ 1 a17 

— z.olDo9/D/9z90oD0bzz90444olzo09 X 10 


9 


+4.09726018688028425780092822833 x 


1 a7 
10' 


1 O 1 OTO OOO tXTO 1Tyl0ylOA0ylA0f?1O1 t^yllO 1 a20 

+Z.I0 / oozoo /ol /4o4o9o49oDlzlo41o X 10 


ID 


—1 .725156234947579334873670891 15 x 


1 a9 
10 


1 O ATAAATI O AOOTtX AA AO C:/^ yl yl A1 1 a23 

-1.09/900/ 1do9oo/o90o9ood4d4oo01 X 10 


ii 


+8.71666539327097376193289896759 x 


1 AlO 

10 


1 1 At^ c;QO0 C^OAi^OTC^O A/l AA/I Q A/l /I At^QO /I A V/ 1 a26 

+i.yooozoozyDz/ooy4yy4oU44Uooz4y X iu 




—5.21094093401193811088710286758 x 


1 a12 
10 


—z.ozi4zyy/yuD/UzyuDii4o/ 00/0090 x iu 




+3.64053240123290947096878491058 x 


1 a14 
10 


1 O ICIAI/lylOyCTI AOO A01 TTi^TOOl AOC^lAv^ 1a32 

+O.lol0144o4/ lyooOzl/ /D/00IO00I9 X 10 


1 A 

14 


—2 .94037039347649323969534446706 x 


1 a16 
10 


A Q c:7A7l CO CQoni 77/5/1 C/l Q/l 7/5n77RQn ^ 1 n35 

—4.00/9/ iDooozUi / /D4o4o4/DO/ /DoU X iU 


15 


+z. / lyo / z4ou / byiiDD/ 4ioUoziyDzoo x 


iU 


1 A COCi A A A TZO OClKn ACiClCl AO CiT^CiA A OaOTZ 1 a38 

+o.4oz044o4oooyo040DD4obo944oDoo X 10 


on 


-1.08008064080683361642422438535 x 


1029 


a 7777700 A o /I o on /I ood o ^5 ccoo 1 a55 
—0. Ill/ MzOddo4oz04dood1oodozood X 10 


30 


-1.60767231445546160409172754060 x 


1053 


— z.oyoool lyo / / / ouoouo / / oyo^^ / ^^4t x lu 


40 


-1.86009993885524622493595668510 x 


1080 


-6.32192703798840441487069765398 x lO^^^ 


50 


-3.13170157718318015959704309179 x 


10109 


-1.86940626135923124849582995810 x lO^^^ 


60 


-2.82752607152079516815460149212 x 


10140 


-2.71176213794760523959296668884 x lO^^^ 


70 


-7.04995099684986774958800182910 x 


10172 


-1.02247285421591824260835479011 x lO^^^ 


75 


+3.31500360451477843327480892738 X 


10189 


+5.81143908209506920438376327293 x lO^^^ 
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Table 1. Continuation A. 



2p — state 


3s — state 


k 


Ek 


77' 


1 


1 ^ ^ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ 

+ 1.50000000000000000000000000000 X 


lO'^ 


+ i.y5^o5i4m^ 111954^ ^zzyyz44oooy X lU 


2 


A r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ r\ 

-4.20000000000000000000000000000 X 


10^ 


—7.99255848804250099004910438109 X 10 


3 


A A f\ f\ f\ f\ r\ e\ e\ r\ r\ f\ f\ f\ f\ f\ f\ r\ e\ 

+4.24000000000000000000000000000 x 


10'' 


+9.86542519182491373199618234956 X 10 


4 


-7.43649066666666666666666666667 x 


10^ 


—2.09268850909736967372029449066 X 10 


5 


+ 1 .87097828977777777777777777778 x 


10^ 


+5.88259296456916749962643192656 x 10 


6 


-6.23210518955750264550264550265 x 


10^" 


A A0/(A1CC1A1A0/1CT'/^01C AAA 1 CT C A^/^ v ^ 1 A 17 

—2.03491551010345668150921550764 x 10 


7 


+2.63615589562176770773494583018 x 


10 


1 O ^ 1 C\ A A O T^AA C ^ CT ^^^H'O C ^^A A O TO O A v > 1 a20 

+8.31244877905745777735770437889 x 10 


8 


-1.38041101245356842308845936964 x 


t r\1 ft 


O A1 AOi? AATTO CT A y< OO/? 1 A1 A y1 OATA^ 1 a24 

—3.91936927735269422610169429796 x 10 


9 


fi t~\ ^^^^ f\f \ f \ ^^^^^ f\ A -1 A tf** f\ A ^^f \ -1 f\ 

+8.78777021937730413143094721050 x 


10^** 


+z.iUoo4zi4lJ7oDzzoD44yioii7oDDOz X iU 


10 


fi \ f\t-\ fy A r~t 1^ f \ ft f\fyf\-t A Ct A f\ -t ft A A fy 

-6.70289485780879092143754218449 x 


-1 r\91 


—1.27856250168620586782771927856 x 10"* 


11 


fi f\ A fy f\ Ct ft -t n ft 1^ JO A A fyCt ft 1^ A J fiftf\f\r\f\ 

+6.04903613875484403654574962209 x 


10^* 


1 O TO A1 A A TO 1 A A ATA A A A A A/^ A C A AT/^ O 1 v ^ 1 a35 

+8.73019078199070904206045207631 x 10 


12 


ft C\ e~\ ft f\ f\ C\ \ f \ ft fi f \ ft ft f\ A fy 1^ fi f \ ^ ft f\f\-t Ck -t ft ft 

-6.38629370869963249560760013166 x 




a i^TTI O C C AA1 /lACAIOOi^AOi^OCOOITCOx^ 1 a39 

—6.67712550914959188602686231758 X 10 


13 


ft f\ ft fv fx f \ ft f \ ^^f \ f% ft A -1 fk f \ H A ft A ^^fk fk 

+7.80769206270237413951487473753 x 




1 C TA1 O /( i^O /I Oi^ /I C OA AA C AAO OOC /I AO C AA 1 a43 

+5.70134624264532995993336408500 X 10 


14 


-1.09524180989414207944239838172 x 


10^4 


C: A 1 A A AATO AA AO O A£;0 O AOT1 £iOf\OGGC\ ^ ^ 1 a47 

—5.41744227892988968803716808669 x 10 


15 


+1. / 4oz/ /4DzoelD4o4ze21oDlzoyUzy X 


iO 


+5.70970876646976596904731808022 X 10 


on 


-9.93004508733968922803863745812 X 


IqSS 


O 11/^O/^CO A A 1AA1 Addl A A CrA11TO/(A1/^v> 1 f\1'2 

—3.11636589919914001005211784216 x 10 


30 


-2.49139318960793926600959408638 X 


1090 


— O.OO ( Oy 100004Zio4iOoUOUDoDlUOOU X iU 


40 


-4.27118446981052432109181315161 x 


10129 


-1.74975503131501618694229226731 x 10^^^ 


50 


-9.94612680182593582345101783873 X 


10170 


-9.32290913905465430025588878862 x lO^^i 


60 


- 1 . 1898340381 1968620569412506832 x 


10214 


-2.21489826289171591863053880143 x 10^62 


70 


-3.81690935673691049129436112647 X 


10258 


-1.28257798832278036045229086676 x lO^^^ 


75 


+2.01879251467438954569448392302 X 


10281 


+2.80538297459908381755301158907 x lO^^o 
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Table 1. Continuation B. 



3p — state 


3d — state 


k 


Ek 




1 


+9.00000000000000000000000000000 X 


lO'^ 


+5.i7i4852o2888045227700755iio05 X 10 


2 


-3.53109375000000000000000000000 x 


-1 r\'^ 

10"* 


— 1. 01742588005743300665089561831 X 10 


3 


+3.44813867578125000000000000000 x 


10^ 


+6.44144210626019957028207894340 X 10 


4 


-5.44958117499880371093750000000 x 


10^ 


—8.64776860916338232383320244317 X 10 


5 


+1.17571125958104512054443359375 x 


lO^'^ 


1 1 C A'\ A^QCCi f\f \ C O AA^C CTIOOA/l^/l/lAIACr v> 1a12 

+1.64145868005800755132947440105 x 10 


6 


-3.18064888985252133830423028128 x 


10^^ 


A 1 r* A 1 AA^cr Or?/^ 1 AO AO AO ^ ^ A C f\c\n <^ A c\ . , 1a15 

—4.16016007587619392285546026342 x 10 


7 


+1.03477309951815151989739404812 x 


10 " 


+1.35506144668665011578479155203 x 10 


8 


-3.95890060949951787171702519397 x 




C A A AAA A COOCiA AT AO 1 OTOTi? AOi^TAA ^ > 1 a22 

—5.45740994533040798137376086700 x 10 


9 


+1.75802363547448200866475465911 x 


10"^' 


1 O O AA /ITOOA/( OOOOOA1 OAT /IT1C010CC 1 a26 

+2.63004738948388301397471581365 X 10 


10 


-8.98508915233072068697058004551 x 


lO-iu 


1 /iTor'ooAi/i/iT/^/iooo/^AAoir'r'A/i/ir'T/i 1 A^iO 

—1.47853891447648836002155944574 x 10 


11 


+5.25232959291187236720839982517 x 


10''* 


1 A COi^OI COCi^OIAOAIA/IAi^l /(OCCOi^OAC 1 a33 


12 


A f\ y^ ft ft y^ e'\ t r\ a ft a ^ r\ ft f^f\ c\ c\ 

-3.49297501706658104776470560233 x 


10""* 


a A C O C O ACCi^ OTAOOC C 1 OT001 OA1 KCil O 1 a37 

—6.95353066637932651872818015613 X 10 


13 


f\ ft f\ f\ A f \ 't ft ft ft fi f \ 't f \ c\ f \ f\ f \ f \ y^ 't ft f \ c\ f \ f \ ft 

+2.62940716668212328895716830967 x 


10*"' 


+5.70612739572397688871373383049 X 10 


14 


't A ft -t ft 't ft ft f\ A ^^f \ f\ ft 't ft ^^f\ fk A A ft 

-2.22927146751916684798616723449 x 


104b 


c: ooc:oAyioooAc:£?£?i 1 yi oto yii£?OA£?y<oc:o 1 a45 

—5.23539432305661143784168964358 x 10 


15 


1 O 1 1 Of\Of\A OO AO AO O A OO A f\C\'l ^1 OaaOOO 

+z.iieUeU4eo4e4eo4oo4UyioioDDzzz x 


1 n50 


+5.34692412323682373417129879084 X 10 


on 


-7.67756482076312202465482219049 X 


10^° 


O CT A A C O OT0 1 AO A/^ A CT A A/^ CT A A O /^O T/^ CT v > 1 A 70 

—2.59058872169806959065242368765 x 10 


30 


-5.13502619416402678854105232212 x 


IQlM 


— Z.Ooo ( 11 ( oyO ( ( DZUDDUlolZUyoDZUU X lU 


40 


-1.78804617789836253655584492158 X 


IQiei 


-1.19414817113404873895159768413 x 10^^^ 


50 


-7.33295913933863525514078947399 X 


10209 


-6.10209745061797542849682305696 x 10^09 


60 


-1.41449923722599152753674938951 x 


10260 


-1.40914942870226252773832277963 x lO^^o 


70 


-6.89053783630996102530832165310 X 


10311 


-7.99404210639204636195381081171 x lO^^i 


75 


+1.39614310998209220877154524815 X 


10338 


+1.73411445429769522609855678396 x lO^^s 



13 



Table 2. Convergence of Pade approximants for energy levels of hydrogen atom in magnetic 
field. The values taken for comparison from Ref. [21] are marked by *\ 



7 


0.1 


0.4 




L 




[L/L-l](7^) 




[L/L-l](7^) 




21 


-0.497526480401260 


-0.497526480401090 


-0.464665 


-0.464599 




22 


-0.497526480401200 


-0.497526480401092 


-0.464658 


-0.464600 


\ls) 


23 


-0.497526480401163 


-0.497526480401093 


-0.464653 


-0.464601 




24 


-0.497526480401140 


-0.497526480401093 


-0.464648 


-0.464601 




25 


-0.497526480401125 


-0.497526480401094 


-0.464644 


-0.464602 


7 


0.16 


0.8 




21 


-0.1236241775347995 


-0.1236241775347930 


-0.098160 


-0.098074 




22 


-0.1236241775347966 


-0.1236241775347925 


-0.098153 


-0.098074 


|2.) 


23 


-0.1236241775347951 


-0.1236241775347928 


-0.098147 


-0.098077 




24 


-0.1236241775347942 


-0.1236241775347928 


-0.098142 


-0.098078 




25 


-0.1236241775347937 


-0.1236241775347928 


-0.098137 


-0.098079 






-0.12362418 *) 


-0.0980892 


7 


0.27 


1.08 




21 


-0.05468786997811 


-0.05468786997805 


-0.045436 


-0.045412 




22 


-0.05468786997796 


-0.05468786997784 


-0.045430 


-0.045398 


\3p) 


23 


-0.05468786997795 


-0.05468786997790 


-0.045428 


-0.045404 




24 


-0.05468786997793 


-0.05468786997788 


-0.045425 


-0.045402 




25 


-0.05468786997791 


-0.05468786997789 


-0.045423 


-0.045404 






-0.05468787 *) 


-0.04540638 *) 



Table 3. Coefficients q of the power corrections to asymptotics of hydrogen atom 

hyper-susceptibilities . 



state 


Cl 


C2 


C3 


C4 


Is 


-2.61829 


+ 1.282 


-2.6 


-11 


2s 


-8.938 


+37.44 


-121 


+2.7 X 10^ 


2p 


-4.6065 


+8.24 


-14.3 


-4 


3p 


-11.227 


+59.5 


-239 


+6 X 10^ 
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FIGURE CAPTIONS 



Fig. 1. Indices of moments of order k that are related by means of the main recurrence relation. 
Each of the links A-D represents one of particular cases. 

Fig. 2. Summation of the PT series for energy levels with the help of Pade approximants. 
£'(7) = [25/25] (7^) - solid curves, -^(7) = [25/24] (7^) - dashed curves. Crosses represent results 
from Ref. [21]. 

Fig. 3. Approach of exact hyper-susceptibilities Ek to their asymptotics for the six states 
of hydrogen atom. 

Fig. 4. Calculation of the ground state wave function at the origin r = by the moment 
method. Solid curve was obtained with the help of Pade approximant [9/9] (7^), dashed curve - 
with the help of [9/8] (7^). Stars denote results of Ref. [23]. 
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